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ABSTRACT 

Context. The detection of the rotational lines of CO in proto-galaxies in the early Universe provides one of the most 
promising ways of probing the fundamental physical properties of a galaxy, such as its size, dynamical mass, gas density, 
and temperature. While such observations are currently limited to the most luminous galaxies, the advent of ALMA 
will change the situation dramatically, resulting in the detection of numerous normal galaxies at high redshifts. 
Aims. We investigate the morphology and strength of the CO rotational line emission emerging from z ~ 3 progenitor 
systems of normal field galaxies seen in the present-day Universe, examine how well ALMA will be able to detect such 
systems and how accurately CO line widths, gas morphologies and ultimately dynamical masses can be inferred. 
Methods. Maps and spectra of rotational CO line emission were calculated from a cosmological N-body/hydrodynamical 
TreeSPH simulation of a z ~ 3 "Lyman break galaxy" of UV luminosity about one order of magnitude below L*. 
To simulate a typical observation of our system with ALMA, we imposed characteristic noise, angular, and spectral 
resolution constraints. 

Results. The CO line properties predicted by our simulation are in good agreement with the two Lyman break systems 
detected in CO to date. We find that while supernovae explosions from the ongoing star formation carve out large 
cavities in the molecular ISM, they do not generate large enough gas outflows to make a substantial imprint on the 
CO line profile. This implies that for most proto-galaxies - except possibly the most extreme cases - stellar feedback 
I effects do not affect CO as a dynamical mass tracer. 

Conclusions. Detecting CO in sub- L* galaxies at z ~ 3 will push ALMA to the limits of its cababilities, and whether a 
source is detected or not may depend critically on its inclination angle. Both these effects (sensitivity and inclination) 
will severely impair the ability of ALMA to infer the gas kinematics and dynamical masses using line observations. 
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?H | 1. Introduction struments (cf. Baker et al. 2004a). With the advent of the 

Atacama Large Millimeter Array (ALMA), this situation 
Since the first discovery of CO in IRAS F10214+4724 at z = may change, making it important that realistic numerical 
2.29 (Brown & Vanden Bout 1992), observations tracking models are dev eloped in order to interpret the data. In spite 
the molecular gas (H 2 ) in distant galaxies have taken on of thiSj the amount of numerical/theoretical work on sim- 
an increasingly important role as a way to study galaxy ulating the expected CO emission from primordial galax- 
formation and evolution (Solomon & Vanden Bout 2005). ies has been limited to a few pi one ering studies (Silk & 
Prime reasons for this are that CO provides an estimate of Spaans 1997; Combes, Maoli & Omont 1999). While these 
the total molecular gas reservoir available for star formation initial efforts were mainly focused on the detectability of 
and black hole accretion and can be used as a tracer of high _ z co lines with future mm-facilities, using 'idealized' 
the dynamical mass of a galaxy (Bryant & Scoville 1996). generic st arburst galaxies as templates, more recently, how- 
The latter is particularly important during the formative everj Narayanan et aL (2006) simulated the CO emission 
stages of a galaxy, where extreme dust-obscuration, as well from tw0 merging gas - r i c h disk galaxies in detail, and stud- 
as outflows, can make dynamical mass estimates derived ied the effects of a massive black hole on the gas dynam i cs 
from optical/near-IR spectroscopy unreliable. of the system as rev ealed by the CO line emission. 

Today, only the most extreme objects at high-z, such as i n this letter we present CO intensity maps and line pro- 

QSOs and luminous submillimeter-selected galaxies (SMGs n i es of a cosmological self-consistent simulation of a merger 

- Blain et al. 2002), are detectable in CO. In contrast, the sys tem, representative of the LBG population at z ~ 3. 

more numerous, but also less massive and presumably less We explore the feedback effects from supernovae (SNe) on 

gas-rich, Lyman break galaxies (LBGs - Giavalisco 2002) the molecular gas, and in particular whether they affect 

are generally too faint to be detected in CO with current in- dynamical mass estimates based on CO. The simulation is 

subjected to the sensitivity and resolution constraints im- 
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posed by ALMA, and we explore the effects this have on 
our results. We adopt a flat cosmology with Q m — 0.27, 
fi A = 0.73 and H = 71 kms" 1 Mpc" 1 (Spergelet al. 2003). 

2. Simulations 

The cosmological simulation of the formation and evo- 
lution of an individual galaxy was performed using the 
N-body/hydrodynamical TreeSPH code, briefly described 
in Sommer-Larsen (2006). The system becomes a Milky 
Way/M31 like disk galaxy at z=0, and is simulated us- 
ing about 2.2 million particles in total, comprising only 
SPH and dark matter (DM) particles at the initial redshift, 
Zi— 39. The masses and gravity softening lengths of SPH 
and star particles are 9.9 x 1O 4 /i _1 M0 and 200/i -1 pc, re- 
spectively. For the DM particles the corresponding values 
are 5.7 x 1O 5 /i _1 M and 370/i _1 pc. The minimum SPH 
smoothing length in the simulation is about 13/i -1 pc. For 
the purposes of this paper, a simulation output at z = 3.26, 
was chosen. We emphasize, that the merger-nature and 
physical properties of the proto-galaxy at this time are de- 
rived self-consistently from a cosmological simulation, in- 
voking detailed baryonic physics, such as radiative gas cool- 
ing, star formation, chemical enrichment and feedback from 
SNe. Note, however, that only atomic radiative cooling is 
invoked, hence neutral, high-density gas does not cool much 
below T ^10 3 K, and "sub-grid" modeling of the molecular 
gas phases are required (see below). Moreover, AGN feed- 
back is not included in the simulations; this may play a role 
in the formation of early type galaxies, but likely not in late 
types (e.g. Sommer-Larsen 2006). 

At z = 3.26, the system has already undergone signifi- 
cant chemical enrichment, and a substantial fraction of the 
gas has assembled into a number of starforming disk sys- 
tems, some of which are in the process of merging. Here, 
we focus our attention on the central 6 kpc of the region, 
which contains the main disk system (A), and a second, 
somewhat smaller system (B) which is undergoing a merger 
with A. Both systems have a combined star formation rate 
of ~ 30 Mq yr" 1 , corresponding to a UV luminosity about 
one magnitude below L* (Sommer-Larsen & Fynbo 2007), 
The stellar masses are M star s — 6.5 x 10 9 M Q (A) and 
M stars — 3.2 x 10 9 Mq (B), which is within the range typ- 
ically observed for LBGs at 2 <, z <, 3 (Shapley et al. 2005). 

Incorporating the formation of molecular gas (H2) from 
the atomic phase (H 1) into numerical models in a fully self- 
consistent manner is currently intractable - largely owing 
to the lack of computational power, and our incomplete 
understanding of how H2 forms in the first place. As a re- 
sult, we adopted physically motivated recipes for the for- 
mation of H2 from the atomic phase. These included the 
expected scaling of the H2 formation rate with the metallic- 
ity (Z) and pressure (P cx n{T + a^)) of the gas (Elmegreen 
1989, 1993), as well as the increase in the H2 destruction 
rate with increasing ambient UV field (G a , in units of an 
average interstellar flux between 6eV <, hv %> 13.6 eV of 
1.6 x 10~ 3 ergem^ 2 s^ 1 - Habing 1968). We ignored, how- 
ever, the destruction of H2 due to cosmic rays (Cazaux 
& Spaans 2004). Also, we made the simplifying assump- 
tion that the formation and destruction of H 2 is in equi- 
librium. In this case, the local molecular gas mass fraction 
can be expressed analytically as a function of P, and 
G (Pelupessy, Papadopoulos & van der Werf 2006), which 
we then used to derive the distribution of H2 in our sim- 



ulation. Within the starforming disks in our simulation, a 
fixed G — 10 5 was adopted, typical of irradiated gas near 
an OB-association. Away from the disks, the average value 
of the local interstellar radiation field was used (G D = 1.7 - 
Draine & Salpeter 1978). The macroscopic velocity disper- 
sion, a v , of the H I gas was calculated from the velocities of 
the SPH particles. We find component A to have molecu- 
lar and atomic gas masses of M(H2) ~ 1.4 x 10 9 Mq and 
M(Hl) ~ 4.0 xlO 8 Mq , respectively, while for component B 
wefindM(H 2 ) ~ 8.0xl0 8 M Q andM(Hi) ~ 1.8xlO 8 M . 
Thus, most of the gas within the disks are in the molec- 
ular phase, consistent with what is seen in local galax- 
ies, where most of the neutral gas resides on much larger 
scales (> 10 kpc). The H2-to-Hi mass ratios are within the 
range observed for luminous starburst galaxies (Mirabel & 
Sanders 1989). In order to resolve scales comparable to that 
of Molecular Clouds Complexes, which in our Galaxy are 
found to have typical size of ~ 50 — 100 pc (Stark & Blitz 
1978), a grid cell size of 10 pc was adopted. Using 5 and 
20 pc cells did not change our findings in any significant 
way. The H2 mass within each cell was assumed to reside in 
a single cloud, following the same large-scale motion within 
the galaxy as the atomic gas, and obeying the same correla- 
tions between velocity dispersion and size/mass as observed 
for a wide range of Galactic clouds (Larson 1981). This then 
allowed us to determine the H2 density as well as the area 
filling factor of each cloud across the grid. The range in 
H2 densities obtained in this way was ~ 50 — 2 x 10 6 cm -3 , 
which closely matches the densities typical of the molecular 
gas responsible for CO line emission in our own Galaxy. 

A Large Velocity Gradient (LVG) model (e.g. Goldreich 
& Kwan 1974) was employed to calculate the CO popula- 
tion levels and source functions for each cloud, using its H2 
density together with fixed values of the kinetic temper- 
ature and CO abundance ([CO/H2] = 10~ 4 ; the C and O 
abundance of the cold gas in the two central disks is approx- 
imately solar) . Guided by studies of local starburst galax- 
ies, which have been shown to harbor a two-phase molec- 
ular ISM (Aalto et al. 1995), clouds with a mean density 
> 10 4 cm~ 3 were given a constant kinetic temperature of 
Tfe = 55 K, while clouds with lower densities were assigned 
a temperature of — 80 K. We used a Cosmic Microwave 
Background temperature of TcMB = (1 + 2) X 2.73 K ~ 
11.63 K, although the effects of the CMB radiation on the 
CO excitation are negligible, except for very cold gas. The 
final emerging CO emission was then obtained by integrat- 
ing the radiative transfer equation along each line of sight, 
taking into account the area filling factor of each cloud. 

3. Results & Discussion 

3.1. Intensity maps and line profiles 

The resulting CO 4-3 intensity maps of the inner 
6kpcx6kpc (corresponding to 0.4" x 0.4" at z = 3.26) 
viewed along the x, y and z lines of sight are shown in 
Fig. [T^-c. The CO 4-3 line is the lowest transition which 
will be accessible to ALMA (band 3) at z = 3.26. It is 
seen that the bulk of the CO emission comes from the two 
merging disk systems, reflecting the fact that the molec- 
ular gas is primarily situated in these two systems, and 
reaches its highest densities there (n(H2) ~ 10 6 cm~ 3 ). 
However, faint CO emission is also seen to originate from 
filamentary structures extending out from the main sys- 
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tern. These are merger-induced tidal streams of diffuse 
(n(H2) ~ 10 1_2 cm. _3 ) molecular gas, hence the low CO 
surface brightness. The diffuse gas accounts for about 1 
and 10 per-cent of the total CO 4-3 and 1-0 luminosities, 
respectively. We find a global CO 4-3/1-0 line luminosity 
ratio of ~ 0.6 from our simulation, which is marginally sub- 
thermal. This, however, is merely an 'aperture effect', due 
to the fact that the extended, diffuse gas contributes more 
strongly to the CO 1-0 than the 4-3 line emission. The CO 
emission shows significant substructure within each of the 
merging galaxies. In the face-on view of the central galaxy 
(Fig. [TJa,), the ISM is seen to consist of multiple cavities of 
low CO emission, surrounded by filamentary regions where 
the CO emission is strong. These structures are created by 
SNe explosions, which carve out 'holes' in the ISM, while at 
the same time creating regions of dense gas as the blastwave 
ploughs through the surrounding medium, compressing it 
in the process. Thus in our simulation, SNe appear to have 
a significant impact on the distribution of molecular gas 
within a galaxy, which is consistent with the notion that 
star formation sites are born in the wake of SNe. 

The CO spectra integrated over the entire region are 
shown in black in Fig. [TJi-f. The first thing to notice is the 
strong dependence of the line shape on the viewing angle. 
As expected, the spectrum in Fig. [TJi is dominated by a 
very narrow (AVpwHM — 50kms _1 ) and smooth single 
profile, which corresponds to the main disk viewed face- 
on. In contrast, the spectra in Fig. QJ-f have double-horn 
profiles, which is normally taken to indicate orbital motion 
in a disk or a merger. While such double-horn profiles ob- 
served towards local ULIRGs have been shown to be due to 
a compact circumnuclear disk of molecular gas (e.g. Downes 
& Solomon 1998), the issue is still controversial for high- 
z objects (Ivison et al. 2001; Downes & Solomon 2003) In 
our case, it is the merger-nature of the system which is 
responsible for the overall double-horn profile. Although, 
we note that each disk system (whose individual spectra 
are shown in red(A) and yellow(B)), when viewed edge-on, 
also give rise to a double-peaked profile, albeit a less distinct 
one. The above illustrates that great care needs to be taken 
when interpreting high-z observations, and underlines the 
need for high-resolution observations capable of resolving 
the emission. Finally, we note that the fact that the line 
profile in Fig. [TJi is smooth, strongly suggest that while the 
molecular gas can be perturbed by SNe, the velocities so 
attained and/or the amount of gas affected are insufficient 
to make a significant imprint on the CO line profile. Of 
course this does not rule out the possibility that in galaxies 
with more prodigal star formation rates, SNe, and possibly 
even stellar winds, could have a significant effect on the ob- 
served CO line profile. Also, we note that Narayanan et al. 
(2006) found that feedback from an AGN in massive merger 
systems can produce outflows of molecular gas, which are 
substantial enough to be detectable in CO. 

3.2. Comparison with observations 

Does the predicted CO properties of the simulated proto- 
galaxy resemble those of observed z ~ 3 LBGs? Only 
two LBGs have been detected in CO to date: the gravi- 
tationally lensed systems MS 1512— cB58 (z = 2.73) and 
J2135-0101 (z = 3.07) (Baker et al. 2004; Coppin et al. 
2007). The two LBGs have star formation rates in the range 
~ 20 — 60 Mq yr , i.e. comparable to our simulated sys- 



tem. The peak CO 3-2 line intensities towards cB58 and 
J2135 are ~ (2mJy)//i and ~ (3mJy)//z, where the lens- 
ing amplification factors are fx ~ 31.8 and ~ 8, respec- 
tively. Thus, their intrinsic CO line fluxes - even if the 
CO 4-3/3-2 brightness temperature ratios differ from unity 
(Papadopoulos et al. 2000) - are within the range predicted 
by our simulation. While the intrinsic CO brightness of 
J2135 is - 7x that of cB58, Fig.QJI-f shows that such differ- 
ences can in principle come about merely from orientation 
effects. Integrating over the entire line, however, we find 
that the total CO luminosity of our simulation is relatively 
constant (within ~ 30 per-cent) with respect to inclina- 
tion angle, suggesting that CO is a robust tracer of the gas 
mass. The CO line widths of cB58 and J213 are very similar 
(AVpwHM — 170 — 190kms _1 ), and in excellent agreement 
with the edge-on CO line widths of components A and B. 
This suggest that both systems are viewed edge-on (J2135 
even has hint of a double-horn profile) and that the differ- 
ence in the CO properties of these two z ~ 3 LBGs is not 
due to inclination effects, but due to an intrinsic diffference 
in the gas content of the two systems. 

3.3. CO as a dynamical mass tracer 

Mapping the CO emission from distant galaxies is one of the 
best ways of determining their dynamical masses, although 
due to the inability of current interferometers to robustly 
resolve the CO emission in high-z objects, most attempts 
at deriving dynamical masses have resulted in little more 
than upper limits. Here we aim to quantify just how well 
one can determine the dynamical mass using CO. 

For each line of sight we have calculated the dynami- 
cal masses of A and B using the disk-formula of Bryant & 
Scoville (1996). Comparing with their total masses (M^ t 
1 .0 x 10 10 M Q and M t B ot = 0.5 x 10 10 M @ ), calculated as the 
sum of the total baryonic mass, i.e. stars and gas (atomic 
and molecular), and DM, we find that in the best case 
scenario when the disks are seen edge-on, the dynamical 
masses can be inferred from the CO emission to within 20 
per-cent. The virial mass of the overall merger-system is 
given by M V i r — 2RV 2 /G, where R and V 2 are the mass- 
weighted distance between the two systems and the line-of- 
sight velocity variance. The actual total mass of the entire 
system (M tot = 2.0 x 10 10 Mq ) is found by integrating up 
all the mass (baryonic and dark) within a sphere centered 
between the two systems and with a radius corresponding to 
half their distance. We find that when the two systems are 
well-separated (Fig. [Tb-c), the total mass can be inferred 
to within 20 per-cent, whereas the discrepancy increases to 
50 per-cent when the two systems overlap as in Fig. [TJi. 

Observations of CO in distant galaxies have to deal 
with limited sensitivity and spatial/spectral resolution, all 
of which can affect dynamical mass estimates. How well can 
we expect to do such observations using ALMA? To answer 
this question we simulated a 24 hr integration with ALMA 
of our system (Fig. Hg-1). We adopted the projected rms 
noise of ALMA at ~ 108 GHz given by the online ALMA 
sensitivity calculator with default settings. The noise was 
added to each 'velocity-slice' in the CO 'data-cube', and 
convolved with a fwhm = 0.1" Gaussian beam - a resolu- 
tion obtainable with ALMA in its extended configuration. 
From Fig. [TJg-l it is seen that only in the case where the sys- 
tem is viewed face-on and is at its brightest, do we see a line 
profile and a detection (4er) in the velocity-averaged maps. 
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Fig. 1. a-c): Velocity-integrated CO 4-3 intensity maps (in units /iJykms -1 ) of our simulation as seen from three 
different lines of sight. Each panel is 6 kpc on a side, with a pixel size of 10 pc. d-f): The corresponding CO spectra 
integrated over the entire field (black curve), the A component (red) and its companion B (yellow), g-i): CO 4-3 maps 
(in units /iJybeam - ) of our simulated system in a-c) after 24hrs of integration with ALMA at 0.1" resolution. The 
maps were obtained by averaging the channel-maps from -200 km s -1 to 200 km s _1 . The contours are — 2, 2, 3, 4 a, where 
a = 10 /xJy beam -1 , where negative contours are shown as dashed curves, j-1): The corresponding spectra (binned to 
20kms~ 1 ) extracted from the centroid position of each source. 



Thus if our simulated system is representative of LBGs at 
z <; 3, then ALMA observations of CO line emission from 
such systems will a) depend heavily on the inclination an- 
gle under which they are viewed, thereby potentially bias 
detections towards face-on systems with relatively small ap- 
parent CO linewidths, and b) essentially be pure detection 
experiments since ALMA will not be sensitive enough to 
constrain CO sizes and line shapes. The implication of a) 
and b) is that it may prove difficult to infer gas kinematics 
and dynamical masses of sub-L* galaxies (such as the one 
simulated in this paper) at high-z using ALMA. 
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